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I.  Summary/Overview 

The  research  objective  of  this  program  is  to  investigate  the  time- varying  dynamics  of  scalar 
mixing  relevant  to  turbulent  reacting  flows.  Specifically,  the  principle  investigator  targeted  the 
development  and  application  of  a  new  class  of  high-repetition-rate  laser-based  diagnostics  to 
measure  the  time-varying  mixture  fraction,  temperature,  species  concentration,  and  scalar 
gradient  fields  in  turbulent  non-reacting  jets  and  non-premixed  flames.  While  the  initial  scope  of 
the  project  focused  on  an  equal  blend  of  measurement  development  and  application,  the  program 
ended  up  heavily  weighted  toward  developing  and  optimizing  new  measurement  capabilities. 
This  was  deemed  necessary  to  ultimately  answer  the  most  fundamental  questions  posed  prior  to 
the  initiation  of  the  project  period. 

The  result  of  the  program  was  five  new  high-speed  imaging  techniques  that  were 
previously  unavailable.  High-speed  planar  Rayleigh  scattering  was  developed  to  measure 
mixture  fraction  and  temperature  fields  in  turbulent  jets  and  flames,  respectively;  high-speed 
spontaneous  Raman  scattering  was  developed  to  measure  the  temporal  development  of  major 
species  concentrations  and  the  mixture  fraction  in  reacting  flows;  high-speed  CH  and  CH2O 
PLIF  were  developed  to  monitor  the  time-varying,  two-dimensional  field  of  combustion 
intermediate  and  radicals  with  varying  levels  of  reactivity. 

Subsequent  to  development,  these  techniques  have  been  used  to  investigate  the  dynamics  of 
turbulent  mixing  and  species  transport  and  will  yield  unique  temporal  and  spatial  statistics  that 
augment  direct  visualization  of  the  scalar  mixing  and  transport  dynamics.  Targeted  applications 
include  understanding  the  time-dependent  interaction  between  the  large-scale  turbulence  and  the 
small-scale  mixing;  the  interaction  between  large-scale  turbulence  and  species 
production/destruction;  and  the  support  of  large-eddy  simulation  (LES)  as  a  method  of  modeling 
time-varying  turbulent  combustion  by  providing  new  temporally-based  visualization  and 
statistical  metrics  for  assessing  model  performance. 

This  work  resulted  in  five  previously-unavailable  imaging  techniques,  seven  archival 
journal  papers,  twelve  conference  papers  and  presentations,  and  eight  additional  invited  seminars 
and  talks.  The  technical  accomplishments  are  expected  to  be  useful  to  the  general  interests  of 
AFOSR  with  specific  interests  expected  in  new  measurement  capabilities  (of  combustion 
dynamics)  that  can  be  applied  to  highly  turbulent  and  high-Reynolds  number  combustion 
systems. 
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II.  Introduction  and  Motivation 

The  United  States  Air  Force  devotes  significant  effort  in  fundamental  propulsion  research  with 
the  goal  of  developing  systems  such  as  scramjets,  pulse  detonation  engines,  and  advanced  gas 
turbine  engines.  The  development  of  state-of-the-art  air-breathing  engines  and  propulsion 
systems  requires  novel  measurement  techniques  capable  of  resolving  governing  physical  and 
chemical  processes  under  highly  turbulent,  high-Reynolds  number,  chemically-reacting  flow 
conditions.  Since  turbulent  combustion  is  the  key  phenomena  in  the  majority  of  Air  Force 
propulsion  systems,  detailed  measurements  under  relevant  conditions  are  critical  for  assisting 
engine/aircraft  designers  in  optimizing  propulsion  systems,  investigating  challenging 
performance  problems,  and  assessing  and  validating  computational  models  such  as  large-eddy 
simulations  (LES)  and  the  underpinning  assumptions  that  could  greatly  reduce  prototyping  time 
and  costs. 

Turbulent  flows  are  inherently  time-varying,  multi-dimensional  phenomena,  and  when 
coupled  to  chemical  reactions,  create  a  highly  dynamic  system  occurring  on  multiple  length  and 
time  scales.  Current  combustion  models  have  difficulty  in  accurately  predicting  the  time- 
dependent  coupling  between  the  turbulent  fluid  mechanics,  scalar  mixing  and  transport,  and 
finite-rate  chemistry  over  a  broad  range  of  conditions.  The  lack  of  robust  and  predictive 
modeling  capabilities  dictates  that  the  primary  assumptions  underpinning  combustion  models 
need  to  be  assessed,  with  carefiilly-designed  experiments  playing  a  crucial  role  in  detailing  the 
combustion  dynamics.  Recognizing  that  turbulence  and  turbulent  combustion  processes  are 
inherently  stochastic,  scalar  mixing  and  transport  processes  must  be  resolved  in  both  space  and 
time  in  order  to  capture  the  dynamic  nature  of  these  flows.  This  requirement  dictates  that  the 
data  acquisition  of  two-  and  three-dimensional  fields  is  at  much  faster  rates  than  typical  time- 
scales  of  the  turbulent  processes  (»  1  kHz).  Non-intrusive  optical  methods  have  the  potential  to 
resolve  pertinent  chemical  and  physical  processes  at  the  smallest  length  and  time  scales;  however 
new  experimental  tools  are  required  to  provide  currently  unavailable  data  that  is  both  spatially- 
and  temporally-resolved. 

Over  the  last  three  years,  the  principle  investigator  has  been  engaged  in  the  development  of 
advanced  high-speed  laser  diagnostic  measurements  in  turbulent  combustion  environments,  with 
a  specific  interest  in  developing  and  applying  new  measurement  techniques  to  examine 
previously-inaccessible  turbulence  and  combustion  dynamics.  As  significant  portion  of  this 
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research  program  was  devoted  to  the  development  of  high-repetition-rate  Rayleigh  scattering, 
Raman  scattering,  and  planar  laser-induced  fluorescence  (PLIF)  diagnostics  to  image  the  time- 
varying  mixture  fraction,  temperature,  and  species  concentrations  fields  in  turbulent  jets  and  non- 
premixed  and  partially-premixed  flames  for  the  first  time.  This  research  laid  the  groundwork  for 
developing  a  better  fundamental  understanding  of  the  physics  governing  the  highly,  unsteady 
turbulent  mixing  and  scalar  transport  processes  within  combustion  environments.  In  addition, 
the  new  measurement  capabilities  support  large-eddy  simulation  (LES)  as  a  method  of  modeling 
time-varying  turbulent  combustion  by  providing  new  structural/dynamic  information  as  well  as 
more  suitable  temporally-based  statistical  metrics  for  assessing  the  predictive  capabilities  of  LES 
models. 

III.  Overview  of  Accomplishments 

The  present  research  program  was  experimental  in  nature,  but  also  resulted  in  synergistic 
activities  with  computational  researchers.  Excellent  progress  was  made  in  the  area  of  high-speed 
imaging  development  and  satisfactory  progress  is  currently  being  made  in  investigating  the 
physical  dynamic  mechanisms  governing  scalar  mixing  and  transport  (in  turbulent  flows  and 
flames)'.  The  research  program  resulted  in  seven  archival  journal  papers,  twelve  conference 
papers  and  presentations,  and  eight  additional  invited  seminars  and  talks.  Highlights  of  the 
previous  accomplishments  are  detailed  in  previous  annual  reports  and  presentations  as  well  as 
the  published  journal  and  conference  papers.  In  this  manner,  only  a  brief  summary  of  the  work 
performed  under  the  present  research  program  will  be  given;  for  a  detailed  description  of  the 
work,  the  reader  is  referred  to  the  individual  published  papers. 

1.  Pulse  Burst-Mode  Imaging 

Recently,  high-speed  imaging  in  combustion  environments  has  been  made  possible  by  advances 
in  diode-pumped,  solid-state  (DPSS)  lasers  and  CMOS-based  camera  systems.  While  the 
emergence  of  commercially-available  laser  systems  has  certainly  been  exciting,  the  low  pulse 
energies  (a  few  milliJoules)  currently  limit  the  diagnostic  techniques  that  may  be  employed  in 
combustion  environments.  Previously,  only  particle  imaging  velocimetry  (PIV)  and  OH  planar 


'  Research  using  the  newly  developed  high-speed  diagnostics  aimed  at  understanding  the  time-dependent  nature  of 
scalar  mixing  and  dissipation  will  continue  beyond  the  current  grant  period  and  at  no  cost  to  the  current  program. 
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laser-induced  fluorescence  (PLIF)  have  been  used  heavily  as  indicated  in  a  recent  review  paper 
by  Bohm  et  al.  [1].  Other  techniques  including  planar  Rayleigh  scattering  or  spontaneous  Raman 
scattering  imaging,  which  are  required  to  measure  important  scalars  such  as  temperature  and 
mixture  fraction  in  turbulent  combustion  environments  [e.g.,  2],  require  very  high  pulse  energies 
and  are  essentially  prohibited  by  the  low  pulse  energies  of  commercially-available  DPSS  lasers. 
Similarly,  many  other  important  combustion  radicals  such  as  CH,  O,  H,  or  NO  cannot  be  easily 
imaged  using  planar  laser-induced  fluorescence  (PLIF)  because  the  UV  laser  energy  obtained 
using  commercial  systems  is  too  low. 

In  order  to  overcome  some  of  these  limitations,  the  PI  and  his  group  use  an  alternative 
approach;  that  is,  pulse-burst  laser  technology  [3]  is  used  to  generate  a  series  of  high-energy  laser 
pulses  in  rapid  succession  to  extend  the  number  of  laser-based  measurements  that  may  be 
performed  at  kHz  (to  MHz)  acquisition  rates.  In  this  section,  a  new  high-energy  pulse  burst  laser 
system  (HEPBLS)  recently  developed  at  Ohio  State  is  described  [4].  It  is  noted  that  this  laser 
system  was  developed  and  constructed  at  no  cost  to  AFOSR,  but  was  leveraged  in  the  PFs 
AFOSR-supported  research  program.  This  system  produces  an  unprecedented  combination  of 
high-repetition-rates,  high  pulse  energies,  and  long  pulse  sequences  to  allow  quantitative  high- 
repetition  rate  imaging  in  turbulent  combustion  environments  that  has  not  been  available 
previously.  This  approach,  which  involves  generating  hundreds  of  sequential  pulses  at  high 
repetition  rates,  is  advantageous  (or  even  necessary)  for  imaging  highly  turbulent  flows. 

A  schematic  of  the  custom  HEPBLS  is  shown  in  Fig.  la.  The  HEPBLS  was  designed  to 
produce  a  unique  combination  of  high  repetition  rates  and  high  pulse  energies  over  long  pulse 
burst  durations  (high  average  power)  for  quantitative  time-resolved  imaging  of  turbulence  and 
combustion  dynamics.  Briefly,  the  HEPBLS  is  a  master  oscillator-power  amplifier  (MOPA) 
system  that  amplifies  the  low-energy  output  (~10  pJ/pulse)  of  a  continuously  running  single¬ 
frequency  Nd:YV04  laser  (which  serves  as  the  oscillator)  in  a  series  of  long-duration,  flashlamp- 
pumped  Nd:YAG  amplifiers.  Between  amplification  stages,  there  are  several  subtle,  but 
important,  system  attributes  aimed  at  controlling  unwanted  birefringence,  diffraction,  and 
thermal  lensing  effects  and  preserving  high  beam  quality,  which  is  important  for  the  high  energy 
and  power  output  of  the  current  system.  Beyond  simple  spatial  filtering,  the  output  of  each 
amplifier  is  “imaged”  onto  the  exit  plane  of  the  subsequent  amplifier  through  relay  optics  and  in 
the  case  of  the  final  stages  (4  and  5),  vacuum  relay  optical  setups.  Each  relay  station  consists  of 
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Fig.  1  -  (a)  Schematic  of  HEPBLS  at  OSU.  (b)  Example  pulse  burst  intensity  trace  and  pulse  energy  distribution 
within  a  single  10-ms  burst.  The  standard  deviation  of  the  energy  fluctuation  is  <  5%. 


a  set  of  lens  pairs  that  focus  the  beam  through  an  additional  “low-pass”  spatial  filter  (diamond 
pinhole  which  is  ~10  times  the  diffraction-limited  spot  size  [5])  and  re-collimates  the  light  at  the 
desired  magnification.  Such  a  process  suppresses  diffraction  ripples  and  helps  to  cancels 
birefringence  effects  at  each  stage. 

After  two  stages  of  amplification,  the  output  energy  (or  more  appropriately  the  intensity)  is 
high  enough  for  efficient  retro-reflection  from  a  Stimulated  Brillouin  Scattering  (SBS)-based 
Phase  Conjugate  Mirror  (PCM)  [6].  The  PCM  serves  to  break  up  what  is  known  as  Amplified 
Spontaneous  Emission  (ASE),  an  undesirable  effect  which  limits  achievable  gain  and  grows 
exponentially  with  the  number  of  amplification  stages.  After  the  PCM  the  pulse  train  is  split  into 
two  parallel  legs  (“A”  and  “B”),  each  of  which  contains  three  additional  dual-flashlamp-pumped 
amplifiers  with  Nd:YAG  rod  diameters  of  9.5  mm  (AMP  3)  and  19  mm  (AMP  4  and  AMP  5). 
AMP3  is  configured  for  double-pass  amplification,  while  AMP4  and  AMP  5  of  each  leg  are 
operated  in  a  single-pass  configuration  with  a  90-degree  quartz  rotator  placed  between  the  two 
amplifiers  to  mitigate  thermally-induced  birefringence.  The  parallel  leg  design  allows  for  a 
delay  in  gain  saturation  in  the  final  amplifier  stages  and  leads  to  a  higher  total  output  energy  for 
equivalent  levels  of  individual  stage  gain.  The  total  gain  for  AMP  3,  4,  and  5  is  greater  than  30 
and  the  total  gain  for  the  system  is  greater  than  10^. 

The  repetition  rate  of  the  master  oscillator  (specified  up  to  250  kHz)  determines  the  pulse 
spacing  within  the  burst,  while  the  fiashlamp  discharge  duration  determines  the  overall  length  of 
the  burst  and  the  number  of  amplified  pulses.  The  current  system  is  designed  to  discharge  for 
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more  than  25  ms  per  burst  and  we  have  tested  the  system  suecessfully  up  to  20  ms.  At  10  kHz 
(100  ps  spacing),  this  corresponds  to  200  high-energy  pulses,  which  is  quite  extraordinary 
considering  that  the  final  energies  are  >  2  Joules/pulse  at  1064  nm  [4]  with  an  RMS  of  the 
pulse-to-pulse  fluctuation  within  the  burst  is  less  than  5%  (see  Fig.  lb).  The  record  lengths  are 
more  than  sufficient  to  characterize  the  time-varying  behavior  of  turbulent  flowfield-flame 
interactions  both  visually  and  statistically  as  well  as  being  well-match  to  typical  LES  simulation 
times  (lO’s  of  milliseconds).  Finally,  the  fundamental,  1064-nm  output  of  the  HEPBLS  is 
frequency-doubled  to  532  nm  for  high-speed  Rayleigh  and  Raman  scattering  imaging  and 
frequency-tripled  to  355  nm  for  high-speed  PLIF  imaging.  Second  harmonic  generation 
conversion  efficiency  as  high  as  60%  has  been  achieved  using  two  type  I  LBO  crystals  (12  x  12 
X  20  mm  )  in  series. 

2.  High-Speed  Mixture  Fraction  Imaging  in  Turbulent  Jets 

The  mixture  fraction  ^  is  the  most  fundamental  (and  perhaps  the  most  important)  scalar  in  non- 
premixed  and  partially-premixed  combustion.  The  turbulent  mixing  field  can  be  uniquely 
characterized  by  and  for  many  combustion  models,  the  turbulent  flame  structure  can  be  related 
back  to  ^  and  its  derivatives.  Practically  speaking,  for  a  large  class  of  combustion  systems,  the 
rate  at  which  the  fuel  and  oxidizer  mix  governs  the  rate  of  chemical  reaction  and  under  highly 
turbulent  conditions,  it  is  highly  desired  to  track  its  evolution  in  both  space  and  time.  A  key 
focus  of  the  experimental  program  was  to  measure  the  mixture  fraction  field,  ^(x,t),  in  turbulent 
non-reacting  jets  to  first  understand  scalar  mixing  dynamics  in  the  absence  of  chemical  reaction 
and  exothermic  heat  release. 

The  mixture  fraction  fields  in  turbulent  non-reacting  jets  were  measured  using  planar 
Rayleigh  scattering.  Rayleigh  scattering  (LRS)  is  a  non-intrusive  diagnostic  that  describes  the 
scattering  from  atoms  and  molecules  whose  effective  diameter  is  much  less  than  the  wavelength 
of  the  incident  light.  For  low-repetition  rate  studies,  planar  Rayleigh  scattering  is  a  widely-used 
imaging  technique  for  mixture  fraction  measurements  in  non-reacting  flows  because  of  the  high 
measurement  sensitivity  that  may  be  achieved  by  issuing  a  fluid  with  a  large  Rayleigh  scattering 
cross  section  into  another  fluid  with  a  much  lower  scattering  cross  section,  i.e.,  propane  into  air 
[e.g.  7-9]. 

The  total  Rayleigh  scattering  intensity  can  be  written  as 
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^RAY  ~  (^) 

where  ^  is  a  constant  describing  collection  volume  and  efficiency  of  the  optical  setup,  Ig  is  the 
incident  laser  intensity,  n  is  the  number  density,  and  is  the  mixture-averaged  differential 

N 

cross-section  defined  as  c  =  V  X  a. .  X  and  o;  are  the  mole  fraction  and  differential  cross- 

mix  i  i  ''  ^ 

i=l 

section  of  species  /,  respectively.  For  a  flow  that  is  isothermal  and  isobaric,  changes  in  the 
scattered  intensity,  Iray  are  due  to  variations  in  amix-  For  a  two-stream  mixing  process  (jet  fluid 
issuing  into  air),  Eq.  (1)  can  be  re-written  as 

i^,=t:iAXj!T,+(i-Xj)<7A  (2) 

where  Xj  is  the  mole  fraction  of  the  jet  fluid  at  any  point  in  the  flow,  is  the  differential  cross- 
section  of  the  jet  fluid,  and  is  the  differential  cross-section  of  air.  Equation  (2)  can  be  re¬ 
arranged  such  that  the  local  mole  fraction  of  the  jet-fluid  {Xj)  is  determined  from  the  Rayleigh 
scattering  signal  {Iray)  at  any  point  in  the  flow  and  measured  calibration  signals  from  pure  jet- 
fluid  (/j)  and  pure  air  (X4): 


^ RAY  ^ A 

’’  I  -I 

T  J-  A 


Subsequently,  the  mixture  fraction  is  determined  from 


X,W,+(l-X,)W, 


(3) 

(4) 


where  7/  is  the  mass  fraction  of  the  jet  fluid  at  any  point  in  the  flow,  and  Yj^  is  the  mass  fraction 
of  the  jet  fluid  originating  from  the  source  {Yj^o=  1),  frj  is  the  molecular  weight  of  the  jet  fluid, 
and  Wa  is  the  molecular  weight  of  the  air. 


Although  frequently  utilized  at  low  repetition  rates,  planar  Rayleigh  scattering  requires 
pulse  energies  that  currently  are  not  available  from  commercial  high-repetition-rate  laser 
systems.  However,  the  new  HEPBLS,  developed  during  the  AFOSR  funding  period,  allows 
high-speed  (multi-kilohertz)  planar  Rayleigh  scattering  imaging. 

For  the  high-speed  mixture  fraction  imaging,  a  turbulent  propane  jet  issuing  into  a  co¬ 


flowing  stream  of  air  was  considered.  The  propane  issued  from  an  8-mm-diameter  circular  tube 
into  a  300  mm  x  300  mm  co-flowing  stream  of  air.  Four  cases  were  considered  where  the 
propane  exits  the  tube  at  Reynolds  numbers  of  10,000,  15,000,  20,000,  and  30,000  based  on  tube 
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Fig.  2  -  25-frame  (of  100),  5-kHz  sequence  of  the  mixture  fraction  field  in  a  tnrbnient  (Re  =  15000)  non¬ 
reacting  propane  jet  issning  into  air.  Image  seqnences  are  shown  for  an  axial  position  of  x/d  =  10.  The 
images  were  acquired  at  10  kHz,  but  only  every  other  image  (200  ps  spacing)  is  shown  for  clarity.  The 
total  time  represented  is  5  ms. 

diameter.  Propane  is  chosen  for  its  large  Rayleigh  scattering  cross-section  (approximately  13 
times  that  of  air),  which  results  in  a  high  signal-to-noise  ratio  for  the  acquired  images. 

As  described  in  Ref  10,  a  considerable  amount  of  effort  was  devoted  by  the  PI  and  his 
group  on  quantifying  the  measured  signals  using  high-framing  rate  CMOS  cameras.  Although 
image  processing  for  traditional  CCD  cameras  is  well-established,  the  characteristics  of  CMOS 
cameras  are  not  well  known  and  this  infrastructure  has  to  be  developed.  In  contrast  to  well- 
characterized  CCD  sensors  which  have  uniform  pixel  response  and  an  exceptionally  high  degree 
of  sensor  linearity,  CMOS  sensors  have  the  potential  to  have  an  independent,  non-linear  response 
for  each  pixel  [11]  because  each  pixel  in  a  CMOS  sensor  acts  as  an  independent  active  circuit. 
Recent  papers  [11,  12]  focusing  on  the  linearity  and  uniformity  of  CMOS  cameras  have  appeared 
in  the  literature  demonstrating  the  need  to  carefully  calibrate  CMOS  cameras  for  quantitative 
measurements.  While  this  may  seem  as  a  minor  point  within  an  experimental  program,  it  cannot 
be  under-estimated  the  effort  that  it  takes  to  characterize  performance  characteristics  of  newer 
technology  such  as  CMOS-based  cameras.  In  this  manner,  the  work  performed  by  the  PI  and  his 
group  in  this  area  will  be  far-reaching,  especially  as  more  and  more  research  groups  transition 
into  high-speed  measurements  and  these  measurements  are  used  to  assess  turbulence  and 
combustion  models  such  as  LES. 

Figure  2  shows  an  example  data  series  taken  from  high-repetition-rate  Rayleigh  scattering 
measurements  of  a  turbulent  propane  jet  (Re  =  15,000)  issuing  into  a  co-flowing  stream  of  air  at 
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Fig.  3  -  Comparison  between  measured  mean  centerline  mixture  fraction,  and  mixture  fraction  scaling 
laws  for  a  turbulent  jet.  (b)  Comparison  between  measured  jet-half-width  (§)  and  scaling  laws  for  a 
turbulent  jet.  Note,  three  Reynolds  number  conditions  are  reported  at  each  axial  position. 

an  axial  position  of  x/d  =10,  where  d  is  the  nozzle  diameter.  The  high-energy  output  of  the  new 
HEPBLS  was  used  at  10  kHz  to  obtain  this  image  sequence.  The  high  signal-to-noise  ratio  (SNR 
>  100  in  the  core)  and  resolved  small-scale  features  are  quite  evident;  in  fact,  the  spatial 
resolution  and  image  SNR  are  sufficiently  high  to  enable  high  quality  determination  of  the  scalar 
dissipation  rate,  %  ~  (V^»VQ,  which  characterizes  the  rate  of  molecular  mixing.  However,  the 
significance  of  the  mixture  fraction  images  lies  beyond  the  spatial  quality,  that  is,  the  utility  of 
the  images  is  evidenced  in  the  temporal  correlation  of  the  images  and  the  ability  to  track 
turbulent  mixing  dynamics  in  real  time.  As  an  example,  the  evolution  of  the  mixing  field 
topology  and  the  dissipation  of  small-scale  scalar  structures  are  clearly  seen,  thereby  displaying 
the  time-varying  nature  of  turbulent  mixing  and  the  steep  concentration  gradients  that  turbulent 
mixing  introduces.  This  work  (as  published  in  [10]  and  an  upcoming  submission  [13]) 
represents  the  first  sets  of  temporally-correlated  (quantitative)  measurements  of  the  mixture 
fraction  fields  in  highly  turbulent  gaseous  jets. 

An  assessment  of  the  accuracy  of  the  mixture  fraction  measurements  is  made  by  comparing  both 
the  mean  value  of  the  centerline  mixture  fraction,  ^c,  and  the  jet  width  for  all  three  Reynolds 
number  cases  and  axial  positions  to  known  scaling  laws  for  turbulent  jets.  The  mean  values  are 
determined  from  200  individual  (non-correlated)  images  (taken  from  the  high-speed  image 
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Fig.  4  -  Example  temporal  correlation  function  for  a  turbulent  (Re  =  15000)  jet  at  x/d  =  15.  The  integral 
time  scale  is  determined  as  -600  ps. 

sequences)  per  Reynolds  number  per  axial  distance.  As  discussed  by  Tacina  and  Dahm  [14],  the 
centerline  value  of  the  mixture  fraction  is  expected  to  scale  as: 

f,(x)  =  5.4(x/rfT'  (5) 

where  x  is  the  axial  distance  downstream  of  the  exit  nozzle  and  d*  is  the  “momentum  diameter” 
which  accounts  for  differences  in  density  between  the  jet  fluid  and  the  ambient  air,  i.e., 
d*={pjet/pair)^^^d,  whcrc  d  is  the  nozzle  diameter.  In  addition,  the  jet-width  (e.g.,  the  “outer 
length  scale”)  can  be  expected  to  scale  as  5(x)  ~  x,  where  the  constant  is  determined  from 
experimental  data  and  the  specific  definition  of  the  jet-width.  Figures  3(a)  and  3(b)  show  that 
both  the  mean  centerline  mixture  fraction  (^c)  and  jet  half-width  measured  using  the  high-speed 
Rayleigh  scattering  diagnostic  compare  very  favorably  to  the  turbulent  jet  scaling  laws  for  all 
three  Reynolds  number  conditions,  thereby  providing  an  initial  validation  of  the  accuracy  of  the 
mixture  fraction  measurements,  at  least  for  the  present  experimental  conditions. 

The  most  recent  work  has  focused  on  using  the  new  high-speed  mixture  fraction  imaging 
capability  to  characterize  turbulent  mixing  dynamics  including  (i)  integral  time  scales  as  a 
function  of  axial  position  and  Reynolds  numbers,  (ii)  power  spectra,  and  (iii)  temporal 
fluctuations  and  correlations  between  mixture  fraction  and  scalar  dissipation  rate.  As  one 
example  of  the  utility  of  these  measurements.  Fig.  4  shows  an  example  of  the  temporal 
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correlation  function  (p)  on  centerline  at  an  axial  position  of  x/d  =  15  for  the  Re  =  15,000  jet.  The 
integral  time  scale  is  calculated  as 

^  =  j  p(T)dT  .  (6) 

which  for  the  present  example  yields  a  value  of  600  ps.  While  this  analysis  can  be  considered 
preliminary,  this  data  represents  the  first  measurements  to  determine  the  integral  time  scale  of 
the  mixture  fraction  field  in  gas-phase  flows  using  optical  diagnostics. 


3.  High-Speed  Temperature  Imaging  in  Turbulent  Flames 

Within  a  turbulent  combustion  environment,  temperature  is  a  very  important  scalar.  For  non- 
premixed  combustion,  the  thermal  field  is  inherently  linked  and  strongly  correlated  with  the 
mixture  fraction  field  in  flame  cases  far  from  extinction.  Thus  thermal  gradients  and  dissipation 
structures  derived  from  the  temperature  field  not  only  yield  physical  information  on  the  thermal 
mixing,  but  on  the  underlying  structure  of  the  molecular  mixing  processes  as  well  [e.g.,  15,  16]. 
In  turbulent  premixed  combustion  systems,  the  temperature  field  yields  a  distinction  between 
reactants  and  products  and  can  be  converted  to  the  characteristic  reaction  progress  variable  (cj), 
which  underpins  many  turbulent  premixed  flame  models  and  yields  direct  insight  into  the 
internal  structure  of  the  premixed  flames.  Furthermore,  local  reaction  rates  and  pollutant 
formation  can  be  strongly  temperature-dependent  and  thus  it  is  important  to  measure  the  local 
temperature  field  and  it  is  highly  desired  to  measure  its  temporal  fluctuation. 

Within  this  program,  the  PI  and  his  research  group  demonstrated  the  first  high-speed 
temperature  imaging  in  turbulent  flames  [17]  and  have  continued  to  improve  the  capabilities 
using  the  new  HEPBLS.  High-speed  planar  Rayleigh  scattering  was  used  to  measure  the  local 
temperature  field  in  the  flames.  For  low  Mach  number  flows,  the  ideal  gas  law,/?  =  nkT,  where 
p  is  the  pressure,  k  is  the  Boltzmann  constant,  and  T  is  the  temperature  can  be  can  be  used  in 
association  with  Eq.  (1)  to  relate  measured  intensity  to  local  gas  temperature: 


^ RAY 


kT 


(7) 


For  flame  temperature  measurements,  the  constant  A  is  typically  accounted  for  by  normalizing 
the  measured  Rayleigh  signal  to  a  Rayleigh  signal  of  known  temperature  and  species 
concentration  such  as  that  from  air  at  room  temperature.  If  the  differential  Rayleigh  scattering 
cross  section  is  constant  across  the  flame,  the  temperature  is  derived  from 
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Fig.  5  -  40-frame  (of  100),  10-kHz  sequence  of  the  temperature  field  in  a  turbulent  (Re  =  22800)  non- 
premixed  jet-flame  issuing  into  air.  Fuel  mixture  is  22.1%  CH4  /33.2%  H2  /  44.7%  N2  and  Re  =  22,800 
(DLR  B  Flame).  Image  sequences  are  shown  for  an  axial  position  of  x/d  =  10.  The  total  time  represented 
is  4  ms. 


T  =  T 


ref 


ref 


'-RAY 


(8) 


where  4e/  is  the  reference  Rayleigh  scattering  signal  from  air  at  room  temperature  {Tref  The 
turbulent  non-premixed  flames  considered  in  this  study  are  the  well-characterized  DLR-A  and 
DLR-B  flames  [18,  19]  which  are  comprised  of  a  fuel  mixture  of  22.1%  CH4,  33.2%  H2,  and 
44.7%  N2.  For  this  fuel  combination,  the  mixture-averaged  differential  cross-section  varies  by 
less  than  3%  throughout  the  flame,  thus  enabling  a  calculation  of  the  local  gas  temperature 
without  a  need  for  measuring  the  local  species  concentrations. 

The  DLR  flames  considered  in  this  study  are  simple  jet-flames  issuing  into  a  low-speed  co¬ 
flow  of  air  (-0.3  m/s)  that  serve  as  benchmark  flames  in  the  International  Workshop  for  the 
Computation  of  Turbulent  Nonpremixed  Flames  (TNF  Workshop)  [20].  The  fuel,  which  consists 
of  22.1%  CH4  /33.2%  H2  /  44.7%  N2,  issues  from  a  8-mm-diameter  tube  at  43.2  m/s  for  DLR 
Flame  A  and  63.2  m/s  for  DLR  Flame  B  into  the  30  cm  x  30  cm  co-flow.  These  flow  conditions 
correspond  to  Reynolds  numbers  of  15200  and  22800,  respectively  based  on  nozzle  diameter. 
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Flame  A  displays  low  levels  of  local  flame  extinction, 
while  Flame  B  is  close  to  blowoff  and  displays 
moderate  levels  of  local  extinction  [18], 

Figure  5  shows  an  example  of  a  recent  10-kHz 
temperature  sequence  obtained  within  a  turbulent  non- 
premixed  flame  operating  at  Re  =  22,800. 

500  mJ/pulse  at  532  nm  was  outputted  from  the 
HEPBLS,  which  allowed  the  use  of  a  CMOS  camera 
only  (no  image  intensifier)  for  increased  signal-to- 
noise  ratio  and  spatial  resolution.  From  the  images. 


Fig.  6  -  Example  temporal  trace  of 
temperature  (at  3  spatial  locatioins)  in 
DLR  flame  B.  Temporal  gradients  as 
high  as  9  X  10®  K/s  are  noted. 


new  features  such  as  the  formation  of  “flame  holes”  (cold  gas  pockets  surrounded  by  hot  gas) 
and  the  upstream  propagation  of  a  high-temperature  reaction  zone  are  clearly  seen,  showing  the 
time- varying  (dynamic)  nature  of  thermal  mixing  and  the  introduction  of  steep  thermal  gradients. 
Image  sequences,  such  as  these,  are  providing  a  new  understanding  of  turbulent  thermal  mixing 
and  flame  dynamics.  Using  Fig.  5  as  an  example,  the  formation  of  the  flame  holes,  the  rapid 
engulfment  by  hot  gases,  and  the  subsequent  upstream  propagation  of  hot  gases  appear  to 
highlight  a  flame  extinction/re-ignition  mechanism.  As  noted  in  Fig.  5,  the  signal  quality 
(signal-to-noise  ratio,  SNR)  is  quite  high  for  these  images.  The  SNR  is  approximately  40  and  17 
for  the  low- temperature  co-flow  and  high- temperature  regions,  respectively.  Figure  6  shows  an 
example  temporal  trace  of  temperature  taken  from  the  data  set  shown  in  Fig.  5  at  three  closely 
spaced  positions  near  the  high-temperature  reaction  zone.  It  is  noted  that  significant  thermal 
fluctuations  are  observed  with  temperatures  fluctuating  as  much  as  900  K  in  100  ps  (9  x  10^ 
K/s). 


To  assess  the  accuracy  of  the  high-speed  temperature  measurements,  the  time-averaged 
temperature  field  and  root-mean  square  (RMS)  of  the  temperature  field  at  an  axial  position  of  x/d 
=  20  are  computed  and  compared  to  the  reference  point-based  data  taken  at  Sandia  National 
Laboratories,  which  appears  as  part  of  the  TNF  database  for  the  DLR  flames  [20].  The  TNF  data 
was  acquired  with  simultaneous  Raman/Rayleigh  scattering  and  CO  laser-induced  fluorescence 
measurements  at  an  acquisition  rate  on  the  order  of  a  few  Hertz  and  spatial  resolution  of  750  pm. 
Approximately  1000  samples  were  acquired  for  each  measurement  location.  For  the  current 
high-speed  temperature  imaging,  the  mean  and  RMS  temperature  values  at  a  given  spatial 
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location  were  extracted  from  the  computed  2D  mean  and  RMS 
temperature  fields  based  on  120  images.  Figure  7  shows  the 
comparison  between  the  current  high-speed  imaging  data  and 
the  reference  data  from  Sandia  for  both  DLR  flames  A  and  B 
plotted  against  the  normalized  radial  position  {r/d).  As  shown 
in  Fig.  7,  there  is  very  good  agreement  between  the  current 
high-speed  imaging  and  Sandia  reference  data  for  the  mean 
temperature  values  for  both  DLR  flames,  both  in  terms  of 
temperature  magnitude  and  radial  profile  shape  (mean 
temperature  gradient).  In  terms  of  the  RMS  temperature 
values,  the  overall  profile  shape  is  computed  quite  accurately 
for  both  flames  A  and  B. 
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Fig.  7  -  Comparisons  of  the  time-averaged  mean  and  RMS  radial  temperature 
profiles  measured  at  an  axial  position  of  x/d  =  20  in  the  DLR  flames  from  current 
high-speed  imaging  (solid  line)  and  the  reference  point  data  (symbols)  from 
Sandia  National  Laboratories  |20].  (Top)  DLR  flame  A;  (Bottom)  DLR  flame  B. 


4.  High-Speed  Raman  Scattering  Imaging  in  Turbulent  Jets/Flames 

A  major  focus  of  the  current  research  program  is  the  development  of  high-energy,  high- 
repetition-rate  diagnostics  that  can  be  used  for  ID  Raman  scattering  line  imaging,  thus 
permitting  a  measurement  of  the  major  species  in  both  space  and  time  and  a  measurement  of 
^(x,t)  in  turbulent  flames  [21].  An  initial  set  of  measurements  demonstrating  the  first  high-speed 
ID  Raman  line  imaging  was  performed  in  turbulent  non-reacting  CH4/H2  jets  issuing  into  air 
[22]  as  shown  in  Fig.  8.  This  flow  was  chosen  such  that  the  major  species  encountered  during 
combustion  were  present  within  the  flow  field.  The  7-mm  line  profiles  of  CH4,  H2,  O2,  and  N2 
are  easily  identified,  although  the  first  measurements  showed  modest  signal-to-noise  ratios. 

Subsequently,  the  PI  and  his  group  have  been  examining  system  performance  within 
combustion  environments.  To  assess  the  accuracy  of  the  high-speed,  CMOS  camera-based 
Raman  imaging  system,  measurements  were  performed  in  a  series  of  well-characterized,  near 
adiabatic  equilibrium,  premixed,  laminar  flames.  Figure  9  shows  the  measurement  of  the  major 
combustion  species  (H2,  N2,  O2,  and  H2O)  and  temperature  over  a  broad  range  of  equivalence 
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Fig.  8  -  10-kHz  sequence  of  ID  Raman  scattering  images  in  a  turbulent  CH4/H2  jet  issuing  into  a  co-flowing 
stream  of  air.  Images  are  centered  at  x/d  =  10. 

ratios.  While  the  measurements  were  performed  at  low  repetition  rates  (10  Hz),  they  were 
performed  with  the  exact  same  optical  and  detection  systems  used  for  10  kHz.  It  is  noted  that  the 
results  match  equilibrium  calculations  very  well  over  the  full  equivalence  ratio  space  giving 
confidence  in  the  ability  to  measure  quantitative  major  species  concentrations  and  mixture 
fraction  profiles  in  turbulent  flames. 

Current  work  in  the  PFs  lab  is  focusing  on  10  kHz  Raman  scattering  measurements  in 
turbulent  H2/N2  flames.  The  significance  of  realizing  multi-kilohertz-rate  imaging  of  ^(x,t)  in 
turbulent  flames  cannot  be  under-stated  as  this  will  give  direct  insight  concerning  the  transient 
mixing  processes  and  dynamic  turbulence-chemistry  interactions  in  reacting  flows  that  directly 
underpin  engine  performance  and  efficiency.  In  addition,  the  new  measurements  of  ^(x,t)  in 
turbulent  flames  will  significantly  augment  the  existing  data  sets  available  to  LES  modelers 
through  temporally-correlated  measurements  and  new  temporal  statistics  such  as  joint 
temperature-mixture  fraction,  temperature -major  species,  major  species-mixture  fraction  space- 


time  correlations. 


Fig.  9  -  (Left)  Photograph  of  a  near-adiabatic  premixed  H2-air  flame.  (Right)  Measurements  (green 
symbols)  of  temperature  and  major  species  as  compared  to  equilibrium  calculations  (red  lines). 
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Fig.  10  —  A  ten-frame  CH  PLIF  image  sequence  in  DLR  Flame  A  at  an  axial  position  of  175  mm 
downstream  of  burner  exit.  The  spacing  between  the  images  is  100  ps,  corresponding  to  a  10-kHz 
acquisition  rate. 

5.  High-Speed  CH  and  CH2O  Planar  Laser-Induced  Fluorescence 

A  final  area  of  research  in  this  program  is  the  application  of  high-repetition  rate  PLIF  imaging  of 
combustion  radicals  and  intermediates  that  cannot  be  performed  using  commercially-available 
laser  systems.  The  output  of  the  PBLS(s)  can  be  frequency-tripled  to  355  nm  for  pumping  an 
optical  parametric  oscillator  (OPO)  as  described  in  Ref  23.  The  development  of  the  PBLS- 
pumped  OPO  has  allowed  the  generation  of  tunable  visible  and  near-IR  laser  light  that  can  be 
frequency-mixed  with  residual  1064  nm,  532  nm,  or  355  nm  output  to  produce  tunable  UV  laser 
light  for  high-speed  PLIF  imaging.  Previously  this  process  has  been  demonstrated  with  NO 
PLIF  imaging  (at  226  nm)  in  non-reacting  supersonic  flows  [24-26]  and  OH  PLIF  imaging  (at 
310  nm)  in  turbulent  hydrogen-air  flames  [27].  In  this  program,  we  have  extended  the 
capabilities  to  cover  CH  PLIF  imaging  at  390  nm  in  turbulent  hydrocarbon  flames  [28].  As  an 
example  of  the  utility  of  this  measurement,  the  CH  radical  has  long  been  considered  an 
appropriate  marker  of  the  inner  reaction  zone  in  both  non-premixed  and  premixed  hydrocarbon 
flames  since  its  location  marks  the  final  step  in  hydrocarbon  fuel  decomposition  [e.g.,  29]. 
Figure  10  shows  an  example  of  a  ten-frame  CH  PLIF  image  sequence  recorded  in  the  DLR 
Flame  A  with  lOO-ps  spacing  between  successive  images,  corresponding  to  a  lO-kHz  acquisition 
rate.  The  images  were  taken  on  centerline  and  at  175  mm  above  the  fuel  tube  with  0.4  mJ/pulse. 
Figure  10  also  shows  two  intensity  profiles  taken  horizontally  from  image  7  of  the  CH  PLIF 
sequence  shown  at  the  axial  location  specified  with  the  two  dashed  lines.  The  SNR  of  the  CH 
PLIF  image  at  both  of  these  axial  locations  is  approximately  6. 

Formaldehyde  (CH2O)  is  an  important  combustion  intermediate  occurring  within  the  “cool 
flame”  region  prior  to  the  primary  reaction  zone.  CH2O  is  a  key  species  in  low-temperature 
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Fig.  11:  10-kHz  sequence  of  the  CH2O  field  in  a  turbulent  non-premixed  CH4  flame  (Re  =  6000).  Images 
are  centered  at  x/d  =  8. 


combustion  chemistry  and  plays  a  crucial  role  in  ignition  and  unwanted  processes  such  as 
“knock”.  Measurements  of  CH2O  fluctuations  allow  an  examination  of  the  effects  of  the  highly 
intermittent  turbulent  flowfield  on  the  low-temperature  chemistry  which  includes  fuel  oxidation 
and  ignition  processes.  An  example  set  of  10-kHz  CH2O  PLIF  images  acquired  using  a  PBLS  in 
a  turbulent  CH4  jet  flame  (Re  =  6000)  are  shown  in  Fig.  1 1  and  were  first  reported  in  Ref  30. 

J 

The  measurements  were  acquired  with  approximately  100  mJ/pulse  at  355  nm  (3  harmonic 
output  of  the  PBLS).  Future  work  includes  using  high-speed  CH2O  PLIF  to  understand  auto¬ 
ignition  dynamics  and  the  temporal  evolution  of  low-temperature  chemistry  in  turbulent  flames 
using  alternative  fuels  such  as  dimethyl  ether. 
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